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Abstract. We study the distribution of the number of (non-backtracking) periodic 
walks on large regular graphs. We propose a formula for the ratio between the variance 
of the number of i-periodic walks and its mean, when the cardinality of the vertex set 
V and the period t approach oo with t/V — > r for any r. This formula is based on the 
conjecture that the spectral statistics of the adjacency eigenvalues is given by Random 
Matrix Theory (RMT) . We provide numerical and theoretical evidence for the validity 
of this conjecture. The key tool used in this study is a trace formula which expresses 
the spectral density of d-regular graphs, in terms of periodic walks. 



1. Introduction 

A graph Q is a set V of vertices connected by a set £ of edges. The number of vertices 
is denoted by V = \V\ and the number of edges is E — \£\. In the present work we deal 
with connected, simple graphs where parallel edges or loops are not allowed. The V x V 
adjacency (connectivity) matrix A is defined such that Aij = 1 if the vertices i,j are 
connected and otherwise. A al- regular graph is a simple graph in which every vertex 
is connected to exactly d vertices. 

Let Qy,d be the ensemble of ci-regular simple graphs on V vertices. Averaging over this 
ensemble will be carried out with uniform probability and will be denoted by E(-). f 
Let P t be the number of t-periodic walks without back-track. It is known that 
E(P t ) = (d—1) 1 for t < \og d _ 1 V, and that the numbers of cycles C t = || are distributed 
as independent Poisson variables [1, 2, 3]. Therefore, 



1 var 




2t , where r = — — > 0. 



(1) 



V' E(P t ) 



f We exclude the set of non-connected graphs since the probability of a non-connected graph is 
exponentially small in Gv,d- 
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We conjecture the following relation valid for any r in the limit t, V — > oo, t/V — > r : 

Fcoe (r) (2) 



1 var(P t ) 



V E(P t ) 

where Fcoe{ t ) is a function derived from Random Matrix Theory for the Circular 

Orthogonal ensemble (COE), and is given explicitly in (31). It takes the following 
asymptotic values: 

( 2r(l + C{d)A + 0(rj) for r -> 

Fcoe(t) = 1 ' . (3) 

[2 for r — > 00 

C(d) = , (d ~ 2) [ - ■ arccoth( v / 2) - — 
V ' ^2d{d - 1) ^ V ; 3tt 

This conjecture stems from our previous work [4, 5]. Here we shall briefly review the 
essence of the previous results, generalize them and present the numerical data which 
substantiate our claim for the validity of (2) and (3). 

To set the scene, we shall define the necessary objects and review some properties 
of (i-regular simple graphs. 



1.1. Definitions 

To any edge b = one can assign an arbitrary direction, resulting in two directed 
edges, e = and e = Thus, the graph can be viewed as V vertices connected 

by edges b — 1, • • • , E or by 2E directed edges e = 1, • • • , 2E (The notation b for edges 
and e for directed edges will be kept throughout). It is convenient to associate with 
each directed edge e = (J, i) its origin o(e) = i and terminus t(e) = j so that e points 
from the vertex % to the vertex j. The edge e' follows e if t(e) = o(e'). 

A walk of length t from the vertex x to the vertex y on the graph is a sequence of 
successively connected vertices x — vi, i>2, • • • , v t — y. Alternatively, it is a sequence of 
t — 1 directed edges ei, • • • , e t _i with o{ei) = Vi, t(ei) = fj+i, o(ei) = £(e t _i) = y. A 
t-periodic walk is a walk of t steps which starts and ends at the same vertex. A walk 
where ej + i 7^ will be called a wa/A; wii/i no back-track or an nb-walk for short. 

In order to count t-periodic nb-walks it is convenient to introduce the 2E x 2E 
Hashimoto matrix Y [6] which describes the connectivity of the graph in terms of its 
directed edges, and avoids back-tracking: 

y e ,e' = <^(e),o(e') — ^e,e' • (4) 

The number of t-periodic nb-walks is P t = try*. 

The spectrum cr(A) = {/J>j}J = i consists of the eigenvalues of A. The largest 
eigenvalue is /ly — d and it is simple for connected graphs. The spectral measure 
(spectral density), from which we exclude the trivial eigenvalue \iy = d, is defined as 

1 V ~ l 
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The spectrum of A is divided in two complementary sets. The first, denoted by R, 
consists of all the eigenvalues which satisfy: \jik\ < 2\/d — 1. The complement, R c 
consists of the eigenvalues for which \fik\ > 2a/ d — 1. The graph is Ramanujan, if 
R c = 0. 

In what follows we shall be interested in the large V limit, and in most cases the 
replacement of V — 1 by V will be justified. We shall do this consistently to simplify 
the notation. 

The trace formula which will be derived in the next section expresses p(/i) as a sum 
of two contributions. The first, often referred to as the 'smooth part' of the spectral 
density, is the celebrated Kesten-McKay measure [7, 8]: 

The second part is called the 'oscillatory part' (or 'fluctuating part') of the spectral 
density. It is an infinite sum over periodic walks on the graph, where each term consists 
of an amplitude which is combinatorial in nature, and some phase. Although this part is 
small compared to the smooth part, it encodes all the interesting features of the graphs, 
including the statistics of the periodic walks. 



2. The Trace Formula 

The starting point for the derivation is the Bass Identity [9] which, for ci-regular graph 
reads: 

det(I (2E) - sY)) = (1 - s 2 ) E - y det(/ (y) (i + (d- l)s 2 ) - sA) . (7) 

The parameter s is an arbitrary real or complex number, I^ 2E "> and are the identity 
matrices in dimensions 2E and V, respectively, and the matrices A and Y were defined 
above. The Bass identity implies that the spectrum of the Hashimoto matrix Y (4) is: 

a(Y) = {(d - 1), 1, +1 x (E - V), -1 x (E - V), 

{Vd^l e^ fc , y/d^l e"^ fe , k = 1, • • • (V - 1))} (8) 

where <f) k = arccos — ^ k , < 9ft(0fc) < n. 
2\/d — 1 

fi k e o-(A) \ {d}. 

We can now write down explicitly trF* which provides the number of t-periodic nb 
walks : 

tr y* = (rf-l)* + 2(rf-l)*/ 2 ( cosh(^fc) + cos(^ fe ) )+i+(£;-y)-(i+(-i)*).(9) 
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where 4>k is defined in (8) and ipk = arccosh , 

It is convenient to introduce the quantities y t , 

ltrY t-(d- iy - 2(d - if' 2 E Mfc6 ^ cosh w k ) 
* = v (dTWi* • (10) 

Not much is known rigorously about the properties of the non-Ramanujan component 
of the spectrum. However, we shall use the known estimates due to Friedman [10] and 
Hoory et. al. [11], to show that y t are bounded as t — > oo. The largest non-Ramanujan 
eigenalue is equal to 2y 'd — 1 • (1 + e), where e is proportional to V~ a and a 0.6. With 
this estimate, the two last terms in the nominator of (10) behave asymptotically as: 

1 



1 + exp ( -t ( ^ \n(d - 1) - C ■ V-°- 3 



(d-iy+2(d-iy /2 J2 cosh(t^ fc ) « (d-iy 

where C is some positive constant. For sufficiently large V 

X - ln(d - 1) - C ■ V-°- 3 > 

and therefore the contribution of the non-Ramanujan eigenvalues is exponentially small 
compared to the leading term, (d — 1)*. Thus, since the leading order term of try* is 
(d — 1)*, the y t are bounded independently of the graph being Ramanujan or not. 
The explicit expressions for the eigenvalues of Y are now used to write, 

-K^) f ^(^) f (-M)')4 ( i:/^) c*, 

where T t (x) = cos (t arccos x) are the Chebyshev polynomials of the first kind of 
order t. An algebraic, straightforward derivation which can be found in [4], results 
in an expression for Pr((j) which is the spectral density restricted to the interval 
M <2 V ^^T) : 

The first term is the smooth part, and can be identified as the Kesten-McKay density. 
We notice that the mean value of yt vanishes as O (^). This is because the counting 
statistics of t-periodic nb-walks with t < log^ V is Poissonian, with E (tr Y l ) = (d— 1)*, 
and because for larger values of t, the leading order term of try* is (d — 1)*. 
As a result: 

Jim E(p(n)) = p KM (n). (14) 

The above can be considered as an independent proof of the Kesten-McKay formula (6). 
The original derivation relied on the fact that (i-regular graphs look locally like trees, for 
which the spectral density is of the form (6). Here, it emerged without directly invoking 
the tree approximation, rather, it appeared as a result of an algebraic manipulation. 



(11) 
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The second term is the aforementioned oscillatory part, which we shall take advantage 
of in the sequel. 

The remainder term in the trace formula (13) is known, and an explicit expression 
is given in [4]. We should mention that (13) is identical to a trace formula derived by 
P. Mnev [12] using a different approach. 



3. Spectral Fluctuations on Graphs and RMT - Introduction and 
Numerical Evidence 

Until recently, the only evidence suggesting a connection between RMT and the spectral 
statistics on graphs was the numerical studies of Jacobson et. al. [13]. In a preliminary 
step in the present research, we performed numerical simulations which extended the 
tests of [13] (see [5]). While describing these studies, we shall introduce a few concepts 
from RMT which will be used in the main body of the paper. 

We shall now (and hereafter), work with the variable <fi rather than /z, as defined in 
(8). This change of variables is well-defined since at this stage we have already taken care 
of all the eigenvalues lying outside the Kesten-McKay support, [— 2y/d — 1, 2y/d— 1]. 
For this reason all the numerics in this paper were carried out using the Ramanujan 
component of c(A), and the number of relevant eigenvalues is modified accordingly by 
V — > V — r c , r c being the cardinality of R c . 

The Kesten-McKay density mapped onto the circle is not uniform: 
2(d - 1) sin 2 

pkmW = - d ^rapy^ • (is) 

The mean spectral counting function provides the average number of eigenvalues up to 
a certain value. It is defined as 

N K m{4>) = V J p K M{(f>)d(t> = ^0 - d -j^- arctan tan ^ . (16) 

Following the standard methods of spectral statistics, one introduces a new variable 
9, which is uniformly distributed on the unit circle. This "unfolding" procedure is 
explicitly given by 

0j = yNKufo) (17) 
The nearest spacing distribution defined as 

P(s) = hm 1e (s - ^{Oj - j , (ft, = Ov), (18) 

is often used to test the agreement with the predictions of RMT (this was also the test 
conducted in [13]). In figure (1) we show numerical simulations obtained by averaging 
over 1000 randomly generated 3-regular graphs on 1000 vertices together with the 
predictions of RMT for the COE [14]. The agreement is quite impressive. 
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S 

Figure 1. Nearest level spacings for 3-regular graphs with 1000 vertices. 
The figure is accompanied with the RMT prediction for the COE. 



Another quantity which is often used for the same purpose is the spectral form- 
factor. This quantity is the main function which we make use of in this paper. It is 
given by 

2 



Kv(t) = 



V 



1 + v E (E cost ( 



0j) 



(19) 



. i<3 



The form-factor is the Fourier transform of the spectral two point correlation function. 
It plays a very important role in understanding the relation between RMT and the 
quantum spectra of classically chaotic systems [14, 15]. 

In RMT the form factor displays scaling. In the limit V, t — > oo ; t/V = r: 

K V (t) = K(T = y) 

The explicit limiting expressions for the COE ensemble is [14]: 

2r - rlog(2r + 1), for r<l 



K C oe{t) 



(20) 



2t+1 



for r > 1 



2-rlog 2r _ 1 , 

The numerical data used to compute the nearest neighbor spacing distribution P(s), 
was used to calculate the form factor, as shown in figure (2). The agreement between the 
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Figure 2. The form factor K(t) (unfolded spectrum) for 3-regular graphs with 1000 
vertices numerical vs. the COE prediction. 



numerical results and the RMT prediction is apparent. This numerical data triggered 
the research in [4, 5]. 

The above comparisons between the predictions of RMT and the spectral statistics 
of the eigenvalues of <i-regular graphs was based on the unfolding of the phases <pj into 
the uniformly distributed phases 9j. As will become clear in the next section, it is more 
natural to study here the fluctuations in the original spectrum and in particular the 
form factor 

/ v 2 

1 



Ky(t) = -E 



5> 



it(f>j 



(21) 



The transformation between the two spectra is effected by (17) which is one-to-one and 
its inverse is defined by: 

6 



<j> = S{9) = N 



KM 



V- 



2tt 



(22) 



This relationship enables us to express Ky(t) in terms of Ky(t). In particular, if Ky{t) 
scales by introducing r = £ then, 

f-vr/2 



V 



71 



(23) 



The derivation of this identity is straightforward, and is given in [5] . 

With this summary of definitions and numerical data we prepared the background 
for the main results of the present work, where we use the trace formula to express 
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the spectral form factor in terms of the variance of the fluctuations in the counting of 
the number of t-periodic nb walks. By assuming that the spectral fluctuations for the 
graphs are given by RMT, we shall derive the variance-to-mean ratio of t-periodic nb- 
walks on graphs. This approach is similar in spirit to the work of Keating and Snaith 
[16] who computed the mean moments of the Riemann £ function on the critical line, 
assuming that the fluctuations of the Riemann zeros follow the predictions of RMT for 
the Circular Unitary Ensemble (CUE). 



4. The Variance-to-Mean Ratio 



The spectral density (expressed in terms of the spectral parameter (8)) is separated 
to its mean and fluctuating parts: 

p R {(t>) = Pkm(4>) + p(4>) ■ (24) 
where Pkm{4>) is defined in (15) and: 

1 

P^) = -y>cos(t0) . (25) 
Using the orthogonality of the cosine, we can extract y t , 

/*7T 

y t = 2 cos (t0)p(0)rf0 (26) 
Jo 

And so: 

E(y*)=4 / cos(t0)cos(#)E(p(0)p(^))rf0# . (27) 
Jo Jo 

From (21), we can write Ky{t) equivalently as: 

K v (t) = 2V / cos(t0)cos(t^)E(p(0)p(^))ci0# , (28) 
Jo Jo 

and comparing (27) and (28) we get: 

K v (t) = (y t 2 ) . (29) 

Since asymptotically, E (P t ) — (d— 1)*, and using (10), equation (29) gives the following 
remarkable equality between a spectral quantity and a combinatorial one, for large V: 

K^(t)-±- E((P t -E(P f )) 2 ) _ 1 V gr(P t ) 
v [ ) ~ 2V' ' E(P t ) ~2V' E(P t ) ■ 1 j 

This is the key ingredient in providing a closed formula for the variance to mean ratio 
of nb-per iodic walks. 

The numerical evidence suggests that the form factor for graphs is given by the COE 
expression (20). Then, combining (30) and (23), we get the desired formula for the 
variance-to-mean ratio of periodic orbits: 

1 var(P t ^ 



v e(p; 



4 / d^pKM^KcoE (-z "~~7T\ ^) = F coe (t) . (31) 

Jo \2irp KM {<p)J 
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The latter is the main result for this paper, providing a formula for the variance-to- 
mean ratio for all values of t. The validity of this result is supported by figure (3) where 
the results of numrical simulations are displayed together with the proposed function 
Fcoe{t). The simulations were carried out by averaging over 100 random choices of 
graphs with V = 1000, d = 10. 
The asymptotic behavior of (31) at the two extremes of r — > and r — > oo, can be 




Figure 3. h 



V< E(p t ) > accom P an i e d by Fcoe(t) (dotted line). 



obtained by using the known behavior of Kcoe- 
At t oo, Kcoe { ^ n ^, ( ^ ) = l - Therefore we get that 

2 



2-Kp KM {4>) 

lim F C oe(t) 



(32) 



This limit is apparent from figure (3). 

At r — > 0, we can expand the middle part of (31) in powers of r (see [5] for details). 
This expansion yields 

Fcoe(t) = 2r ■ (1 + C(d)V^ +■■■) (33) 
where C(d) is given explicitly by (3). 

The most striking feature lies in the fact that the deviation from the Poissonian 
expression is of order r 1//2 . This is illustrated in figure (4) where ■ — 2r) / (2C(d)) 
is plotted for graphs with various values of d. The expected power-law and data collapse 
are clearly visible for r < 0.2. 



In our opinion, the numerical evidence presented above is convincing enough to 
suggest that our main conjecture (2) is valid. A combinatorial approach is called for to 
test it further. 
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